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Abstract 

The effect of Rician-ness on the capacity of multiple antenna systems is investigated under the assumption that 
channel state information (CSI) is available only at the receiver. The average-power-constrained capacity of such systems 
1 is considered under two different assumptions on the knowledge about the fading available at the transmitter: the case 
' in which the transmitter has no knowledge of fading at all, and the case in which the transmitter has knowledge of 
the distribution of the fading process but not the instantaneous CSI. The exact capacity is given for the former case 
' while capacity bounds are derived for the latter case. A new signalling scheme is also proposed for the latter case and 
' it is shown that by exploiting the knowledge of Rician-ness at the transmitter via this signalling scheme, significant 
capacity gain can be achieved. The derived capacity bounds are evaluated explicitly to provide numerical results in 
' some representative situations. 
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I. Introduction 



Prompted by recent results suggesting possible extraordinary capacity gains P, multi- 
'ple transmit and receive antenna systems have received considerable attention as a means of pro- 
^ viding substantial performance improvement in wireless communication systems. In such multiple- 
y—^ ! input/multiple-output (MIMO) systems, multiple transmit /receiver antenna combinations provide 
^ 'spatial diversity by exploiting channel fading. It has been shown in [2], [3] that when the receiver 
I has access to perfect channel state information but not the transmitter, the capacity of a Rayleigh dis- 
O . tributed flat fading channel will increase almost linearly with the minimum of the number of transmit 
^ [ and receive antennas. 

' In most previous research on the capacity of multiple antenna systems, however, the channel fading is 
O assumed to be Rayleigh distributed. Of course, the Rayleigh fading model is known to be a reasonable 
J> . assumption for fading encountered in many wireless communications systems. However, it is also 
^ ! of interest to investigate the capacity of multiple antenna systems when the Rayleigh fading model 
5-H is replaced by the more general Rician model. Not only does this generalize the previously derived 
■ — ■ capacity results, since both additive white Gaussian noise (AWGN) and Rayleigh fading channels may 
be considered to be limiting cases of the Rician channel, but Rician fading is also known to be a better 
model for wireless environments with a strong direct Line-Of-Sight (LOS) path In this paper, we 
consider the capacity of multiple antenna systems in Rician fading for two cases of interest: that in 
which the receiver has perfect channel state information (CSI) but the transmitter has no knowledge 
of the fading statistics; and that in which the receiver has perfect CSI and the transmitter knows the 
distribution of the fading process, but not exact CSI. 

We begin, in Section |nj by introducing the multiple antenna system model of interest and the 
assumptions on the fading process. Next, in Section UTTl we address the general capacity problem for 
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the Rician fading channel and obtain an upper bound for the capacity of muhiple antenna systems 
under Rician fading. We exphcitly evaluate this upper bound for some special cases. 

In Section IIVI we investigate the exact capacity of MIMO systems in Rician channels under the 
assumptions of perfect channel state information at the receiver and no knowledge of the fading distri- 
bution at the transmitter. In addition to providing a lower bound on the capacity for the case in which 
the transmitter does know the fading distribution, this result also serves as a measure of the capacity 
variation of a system designed under Rayleigh fading assumption but operating in an environment 
where a strong LOS component is present. This is because the capacity-achieving distribution for 
the case considered in this section also achieves the capacity in the Rayleigh channel. We explicitly 
evaluate this capacity for some interesting special cases. 

We will see that there is a large capacity gap between the upper bound of Section lllll and the lower 
bound of Section |TV| obtained with signals designed to be optimal for Rayleigh fading. In Section |V| 
we propose a new signalling scheme for multiple antenna systems with perfect-CSI at the receiver and 
only knowledge of the fading distribution, but not the exact CSI, at the transmitter. We derive tight 
upper and lower bounds for the capacity of a multiple transmit antenna system with this new input 
signal choice. By comparing these capacity bounds with the results obtained in Section IIVI we will 
show that the exploitation of the knowledge of the fading distribution at the transmitter can provide 
significant capacity gains. Finally, we finish with some concluding remarks in Section IVTl 

Some mathematical results that we will need in the rest of the paper are given in the Appendix. 

II. Model Description 

We consider a single user, narrowband, MIMO communication link in which the transmitter and 
receiver are equipped with Nt and Nn antennas, respectively. We consider the ideal case in which 
the antenna elements at both transmitter and receiver are sufficiently far apart so that the fading 
corresponding to different antenna elements is uncorrelated. The discrete-time received signal in such 
a system can be written in matrix form as 

y(z) = H(^)x(^) + n(0 , (1) 

where y{i), x(i) and n{i) are the complex iVij- vector of received signals on the receive antennas, the 
(possibly) complex NT-vector of transmitted signals on the Nt transmit antennas, and the complex 
iVfi-vector of additive receiver noise, respectively, at symbol time i. The components of n{i) are 
independent, zero-mean, circularly symmetric complex Gaussian random variables with independent 
real and imaginary parts having equal variances; i.e. n{i) ~ ^^0(0,!^?^), where Itv^ denotes the 
Nr X Nji identity matrix. The noise is also assumed to be independent with respect to the time index. 

The matrix in the model (Q) is the Nji x Nt matrix of complex fading coefficients. The (n^, n^)- 
th element of the matrix denoted by (H(i))^^^^, represents the fading coefficient value at time 

i between the riR-th receiver antenna and the n^-th transmitter antenna. The fading coefficients in 
each channel use are considered to be independent from those of other channel uses, i.e. {H(z)} is an 
independent sequence. As noted in , this gives rise to a memoryless channel, and thus the capacity 
of the channel can be computed as the maximum mutual information, 

CNT,Na = max X(x;y) 

where Px is the probability distribution of the input signal vector x that satisfies a given power 
constraint at the transmitter and X(x; y) is the mutual information between the input x and output 
y. In this case, we may also drop the explicit time index, i, in order to simplify notation. 

The main purpose of this paper is to extend the previously known capacity results for multiple 
antenna systems in Rayleigh fading to Rician channels. Thus, we will assume that the elements of H 
are Gaussian with independent real and imaginary parts each distributed as J\f{fi/^/2,a'^). Moreover, 
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the elements of H are assumed to be independent of each other. So, the elements (H)^^ of H 
are independent and identically distributed (i.i.d.) complex Gaussian random variables (H)^^ ~ 

Mc (^"^(l + i)) 2(7^ j , for = 1, . . . , Nr and = 1, • • • , iVr, and the distribution of the magnitudes 
of the elements of H have the following Rician probability density function (pdf): 

fair) = 2(1 + ^)re-(^+'^)'-'-%(2v/«:(l + K)r), (2) 

where /q is the zero'th order modified Bessel function of the first kind j^J and we have introduced the 
Rician factor, k, defined as 

^ = 2^- 

For notational convenience, we have also introduced the normalization + 2cr^ = 1. Note that (j21) 
reduces to the Rayleigh pdf when k = (which implies that /i = 0). 

When elements of H are distributed as described above we say that H is a complex normally 
distributed matrix, denoted as H ~ Mc (M,IiVj, ® S) where S is the Hermitian covariance matrix of 
the columns (assumed to be the same for all columns) of H and M = i^jH}. For the assumed model. 



and 



2^'Iiv«, (4) 



M = ^(l+j)*iv«,7v^, (5) 

where ^'at^^atj, denotes the Nr x Nt matrix of all ones. 

Next, let us define n = meix{NR, Nt}, m = mm{Nf{, Nt} and 

/ HH^ if NrKNt 
W = I H^H ifNn>NT ■ 

Then, W is always an m x m square matrix. It is known that when H is a complex normally 
distributed matrix as described above, the distribution of W is given by the non-central Wishart 
distribution [HI, [ZI, |H] with pdf 

/w(W) = e-*H^"^^^"}oFi(n;S-iMM^S-^W)/°,(W), (7) 
where in (jTj) /w(W) denotes the (central) Wishart pdf: 



/O (W) = g-trS-iW|^|n-m .gN 



which results when the elements of H are iid zero mean Gaussian random variables, and where the 
complex multivariate gamma function and the Bessel function of matrix argument o-^i are defined 
in the Appendix (see ^) and ^^ V. 

Note that in (jZj) we have assumed, without loss of generality, that W = HH-'^. We will continue 
to use this assumption throughout unless stated otherwise. We use the shorthand notations W ~ 
[n, S) and W ~ [n, S, S^^MM^) to denote that W has the Wishart distribution with pdf 
(jHl) and that W has the non-central Wishart distribution with pdf (|7j), respectively. 

^Note that (gSJ applies in this case by noting that S^^MM^S^^W = E^^MM^S^^HH^ and the trace relationship 
tr(AB) = tr(BA). 
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III. Capacity of the Multiple Antenna Rician Fading Channel 

It is shown in j2] that, for the Rayleigh flat fading channel (i.e. the model of Section |n] with k = 0) 
under the total average power constraint E{x.^x.} < P, the capacity of the channel is achieved when 
X has a circularly symmetric complex Gaussian distribution with zero-mean and covariance -^Inj,, 

and that this capacity is given by the expression C^^^^^ = i?|logdet (Inr + i^HH-^j |. 

However, it is also easily shown that the capacity achieving transmit signal distribution x for a 
multiple antenna system under the average power constraint above is circularly symmetric, zero-mean 
complex Gaussian regardless of the actual fading distribution, as long as the receiver, but not the 
transmitter, knows the channel fading coefficients. Thus, only the covariance matrix Q = E{x.:k^} of 
the capacity-achieving distribution depends on the fading distribution, and the capacity of the multiple 
antenna system is then given by 

CNr,Nn = max EH{logdet(HQH^ + I^^)} . (9) 

tr Q < P 
q > 

In the case of deterministic fading where the matrix H has all its elements equal to unity (i.e. the 
Rician model with k — > oo) and this is known to the transmitter, the so called water-filling algorithm 
jH] specifies the covariance matrix structure of this capacity achieving Gaussian distribution to be of 
the form Q°° = ■^'^^^ where denotes the Nt x Nt matrix of all ones 0. In this case, the 
capacity is given by 

C^r,N, = \ogil + NnNTP). (10) 

Alternatively, in the case where the fading is Rician but without knowledge of k at the transmitter, 
the capacity achieving distribution is the same as in the Rayleigh case P, [2], [Sj, i-e. its covariance 
matrix is Q° = -£^1^^- 

Thus, for a channel with Rician distributed fading having a general value of k, which is known to 
the transmitter, one would expect the covariance matrix Q of the capacity achieving distribution to 
lie in between these two extremes. Although the capacity-achieving Q for this case is unknown, in the 
following paragraphs we derive an upper bound for the capacity of this channel. 

Observe that for any Q the matrix HQH^ + I^^ is positive definite, and that the function logdet 
is concave on the set of positive definite matrices. Thus, applying Jensen's inequality to we have 

Cnt,Nr < max logdet (iV^QT + Iat^) , (11) 

tr Q < P 
q > 

where we have used the determinant identity det{I + AB) = det{I + BA) and introduced the notation 

^{H^H} = NrT. (12) 
It is easy to show that the x A^^ matrix T is given by 



1 + K 



1 + K K ... K 

K, 1 + K, . . . K 

K, K, . . . 1 + K 



(13) 



We observe that for any k, such that < k < oo the matrix T is non-singular and thus all the 
eigenvalues of T are non-zero. In fact, if we denote the eigenvalues of T by Aj for i = 1, . . . , Nt, then 
it can be shown that, 



1+Ntk \r j = 1 
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We may decompose Y as T = UDU'^ where D is the Nt x Nt diagonal matrix having the 
eigenvalues in (jl4|) as its diagonal entries, and U is a unitary matrix. Substituting this into the right 
hand side of (fTT|) . we have 



Cnt,Nr < max ^ log det ( NrQD + Int j , 

q >~0 



(15) 



where we have let U'^QU = Q. Now it is easy to see that the right-hand side of (fT3j) is maximized 
by a diagonal Q and the diagonal entries are again given by the well-known water-filling algorithm. 
Indeed, one can show that the maximizing diagonal matrix Q is such that. 



mm 



P k{1+k) 

P _ 

Nt 



Nt + 


P 

Nt 


k{1+k) 





k(1+k) 
Nii(l+NTK) 



Nr{1+Ntk) 



if i = 1 
Hi = 2 Nt 



for < K, < oo. 



where [x]^ = max{0,x}. 

Thus, for < K < oo, the capacity of a multiple antenna system in Rician fading, subjected to the 
average transmit power constraint P, with perfect CSI at the receiver and knowledge only of k at the 
transmitter is upper bounded as 



C 



Nt,Nii 



< log 1 



Nn (1 + Ntk) 



mm 



P 



k(1 + k) 



+ (Nt - 1) log 1 + 



N 



R 



Nt' Nr{1 + Ntk) 
P k(1 + k) 



Nt + 



P 

N^ 



k(1 + k) 



Nn{l + NTK) 



1 + K 



Nt Nr{1 + Ntk) 



(16) 



(For K — > oo, the exact capacity is given by (fTTH) .) 
Next, we will illustrate the above bound for some special cases. 

Case 1: k = 

It is easily seen that for k = 0, which is the Rayleigh fading case, the above bound reduces to 

Nr 



c 



Nt,N, 



^ < iVrlog 1 



-p 



for K. 



0. 



(17) 



From ()17|). we see that when Nr = Nt the capacity upper bound is a linear function of Nt- In 
fact, it was shown in j3] that in this case the capacity can be approximated by a linear function of Nt 
asymptotically for large numbers of antennas. 

In addition, if Nt = 1 in ()17j) . then 



Cnt,n^ < log(l + NrP) 



for K = and Nt = 1. 



In fact, it was shown in that the capacity of the Rayleigh fading channel in this case is asymptotic 
to log(l + NrP) for large Nr. 



Case 2: Nt = 1 
When Nt = 1, the capacity upper bound in (jl6|) becomes 

Ci,Nn < \og{l + NrP) for K > 0. 



(19) 



From (jl9|) we see that the bound (jl8|) is in fact valid not only for the Rayleigh channel but for any 
Rician channel with an arbitrary value of n, in this case of Nt = 1. 
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Case 3: Nr = 1 
The capacity of the Rician channel is bounded in this CclSG clS 

P K(l + h 



Cnt,i < log I 1 + ^ _^ ^ I mm 



+ (Nt - 1) log 1 



Nt 1 + Ntk 



Nt + 



P + 



1 + Ntk 



1 



1 + K 



P k{1 + k) 
Nt 1 + Ntk 



(20) 



Case 4: Nr = Nt = n 
In this case, the capacity upper bound of (fTT)|l reduces to 

Cn,n < log 1 + — I min <! P, ^ >- } n + 



K 



+ (n - 1) log 1 + 



1 + riK 



P 



P 



k{1 + k) 
1 + nn 



K 



k{1 + k) 
1 + nn 



Note that the exact capacity as k, — > oo equals log(l + n'^P) in this case. 

IV. Capacity of the Rician Channel When the Transmitter Does Not Know The 

Fading Distribution 

In this section, we investigate the capacity of the average-power-constrained Rician channel when 
the receiver (but not the transmitter) has perfect CSI, and the transmitter does not know the fading 
distribution. Of course, this capacity provides a lower bound for the capacity in the situation where 
the transmitter does know the fading distribution. Recall from the preceding section that the opti- 
mal transmitted signal distribution when the transmitter does not know the fading distribution is a 
circularly symmetric complex Gaussian distribution with i?{xx^} = -^Ijvy. Since this distribution is 
also optimal for the average-power-constrained Rayleigh channel, evaluation of the capacity with this 
signal distribution in the Rician channel also quantifies the capacity variation of a system designed to 
be optimal for Rayleigh fading but operating in a Rician channel (i.e. a channel with a line-of-sight 
component). 

Applying this signal distribution, from (jH)) the capacity of the Rician channel (arbitrary k) with no 
transmitter knowledge of n is given by 



C 



Nt,Nii 



Ew<ilogdet ( ^W + I„ 



(21) 



where the expectation is with respect to W ~ Wm {n, S, S ^MM-'^) with pdf (jlj). 

Since W in ()2H1 is an m x m Hermitian matrix, if we denote its (non-negative) eigenvalues by 

Ai, A2, . . . , Am, then we have det ^-^W + j = YYiLi (^^ + t^-^^) • Hence, the capacity in (f?T|l can 
be given in terms of the eigenvalue distribution of the non-central Wishart distributed matrix W as 



E 



Ai,...,A„ 



logHfi 



i=l 



mEx <! log ( 1 + ^A 




(22) 



(23) 



where in (I23p we have taken A to be any un-ordered eigenvalue of the non-central Wishart distributed 
random matrix W. 
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A. Special Case 1: Capacity in the Limit of Large 

Before attempting to evaluate the capacity exactly, it is instructive to investigate its behavior in 
the limit as the number of transmit antennas increases without bound. This will also allow us to 
compare the asymptotic capacity in this Rician case with that of the Rayleigh case given in |3j, 
thereby illustrating the effect of the non-zero mean of fading coefficients on the capacity. 

Note that for a fixed A^^^;, the elements of the matrix HH'^ are the sums of iid random variables 
with finite moments, and thus by the strong law of large numbers (SLLN) we have, almost surely, 

T where the matrix T is defined in (jl3|) (it is now taken to be an m x m square 



lim 



Nt- 



matrix). Thus, for a fixed number of receive antennas, when the number of transmit antennas becomes 
very large, the capacity of this channel is given by 



a 



lim logdet (Itv^ + ^HH 

Nt — ►oo \ J\rp 



H 



(Nn-l) log 



1 + 



P 



1 + K 



+ log 



logdet (I^^ + PT), 
P 



1 + {Nnf^ + 1) 



1 + K 



(24) 



The following two cases illustrate the dependence of the above asymptotic capacity on the Rician 
factor k: 



lim Coo,iVH = iV^log(l + P) (25) 

K >0 

and 

lim Coo,iv« = log(l + iV^P). (26) 

K >00 

Note that is the asymptotic capacity of the Rayleigh channel given in It is easily seen that 
Coo,Nji in is a monotonically decreasing function of k for k > and > 1 and is constant for 
= 1. Thus, we observe that when Nt is very large and N^ > 1, increasing the determinism of the 
channel lowers the capacity if the transmitter is not aware of this increase, and moreover the Rician 
fading environment will degrade the capacity of a system that is designed to achieve the Rayleigh 
channel capacity. Of course, this does not necessarily mean that the capacity of the Rician fading 
channel is less than that of the Rayleigh fading channel when the transmitter knows k. In fact, from 
fllUp we may recall that in the deterministic case, which is the limiting case of Rician fading with 
K — > oo, the capacity of the channel, as achieved by the water-filling algorithm, is known to be 
log(l + NrNtP), for any N^ and Nt- Similarly, it is reasonable to expect that the capacity of a 
multiple antenna system in Rician fading, with an arbitrary n, to be greater than that in a Rayleigh 
fading channel when transmitter knows the value of k. 

Figures Q and |21 show the dependence of the asymptotic capacity (j24j) on the Rician factor k for 
P = dB and P = 10 dB, respectively. These graphs illustrate our conclusion about the asymptotic 
capacity degradation of the Rician channel with increasing k when the transmitter has no knowledge 
of K. 

In Fig. 121 we show the asymptotic capacity versus the number of receiver antennas Nn for 
K = 10. Figure El shows the almost linear dependence of this capacity on the number of receiver 
antennas (which is the smaller of Nt and Npt in this case), similarly to the previously established 
linear dependence for the Rayleigh fading environment (PP, [SI). 

B. Special Case 2 .■ min{Nji, Nt} = 1 

In this section we will evaluate the exact capacity of the Rician fading channel (subject to the signal 
choice assumed throughout this section) in the special case m = mm{Nfi, Nt} = 1. In this special 
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case (IH) reduces to a scalar S = 2cr^, and thus the pdf of W in ((Tj) (which is a scalar) can be written 
as 



/wW = e-^"-^^^e-^|iy|"-i + k)w) ^" '^/„_i(2v/n/€(l + /t)iy), (27) 

where we have used the fact that Ti{a) = T{a) and (j46p from the Appendix. 
From the capacity in this special case is 

CNr,Nn = j^" log (^1 + MW)dW, (28) 

where fw{W) is given in ()27p above. As noted by Telatar in jS] for the Rayleigh fading channel, from 
()28p we observe that the capacity is not symmetric in and Nt also in the Rician case. Thus, we 
have two cases to consider as below. 

B.l Case 1: Nr>Nt = 1 

From (j27p and (|28p. the capacity of the Rician fading channel when the transmitter does not know 
K in this case is, 

1 f°° 

Ci,;v« = ^pTTTT / log{l + PW)W''--'e-'^^{W,Nn)dW, (29) 

where we have introduced the function 

m,n) = ^ ^^"^ ^, (^)"e-^iy-^^"-^^4-i (2Vn.{l + .)w) . (30) 
[nK(l + fi;)]2(""^) \ J ^ ^ 

Figures 0] and El show the capacity of the Rician fading channel when the transmitter does not know 
K in this special case of = 1, against the number of receiver antennas for k = 1 and k = 10, 
respectively. Included on the same graphs are the corresponding capacity curves for the Rayleigh 
fading channel (k = 0). We observe that the capacity of the Rician channel is greater than that of 
the Rayleigh channel and the capacity gap increases with increasing values of k. We can also observe 
from Figs. |3] and El that the capacity gap is prominent for smaller numbers of receiver antennas and, 
as — > oo, the two capacities converge to the same value. In Figs. 01 and El we have also shown the 
capacity upper bound for this system given by (fT^. It is clear for these figures that the upper bound 
f)19p is very tight in this case. Indeed, since Nt = 1, in this case the optimal signal covariance that 
satisfies the average power constraint is Q = Q = E{x'^} = P. 

B.2 Case 2: A^'t > A^^ = 1 

The capacity of the Rician fading channel when the transmitter does not know k in this case is 



(31) 



where ip{W, Nt) is given by (jHUj) . 

Figures El and [71 plot the capacity of the Rician fading channel in this special case of A^r = 1, against 
the number of transmit antennas for k = 1 and k = 10, respectively. As before, included on the same 
graphs are the corresponding capacity curves for the Rayleigh fading channel (i.e. the k = case). 
Again, we observe that the capacity of the Rician channel is greater than the capacity of the Rayleigh 
channel and the capacity gap increases with increasing values of k before finally converging to the same 
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value for large A^^. Particularly, from Fig. [7|we note that for a smaller number of transmit antennas 
the capacity gap is significant. 

Shown also on Figs. El and [7| is the capacity upper bound for the A''^^ = 1 case given by (pUjl . 
From these figures we observe that, unlike in the case of Nt = 1, the upper bound is very loose for 
the case of = 1. However, recall that we are using a particular input signal distribution which 
is not necessarily the capacity achieving distribution for this particular channel under the assumed 
conditions on CSI and fading statistics. Rather, we were assuming a signal distribution that is only 
optimal for the Rayleigh fading channel or for a system in which transmitter does not know the value 
of K. Thus Figs, ini and [7| suggest that scaled identity matrix might not be the form of the covariance 
matrix of the capacity achieving input signal distribution for a multiple antenna Rician channel when 
the transmitter knows k, and with better signal choices that exploit the Rician-ness inherent in the 
fading distribution, we may be able to obtain higher capacities. In Section |V] we propose a new choice 
for the covariance matrix which offers much higher capacity than that achieved by the scaled identity 
matrix. 

Using a series representation of the modified Bessel function, it is straightforward to show that 

lim^ ,Qilj(W,m) = 1. Thus, in the limit n — > 0, (PIJ|) and reduce to the corresponding capacity 

expressions for the Rayleigh channel given in [3^, as one would have expected. 

C. General Capacity Expression for the Rician Channel 

In order to compute the capacity of the Rician channel for an arbitrary number of transmit / receiver 
antennas, as given in (|23p. we need to find the latent root distribution of the non-central Wishart 
distributed matrix W. This latent root distribution has been studied previously (|5j, [TU], [TT], [7j) 
and, in particular, we have the following result from [7j. 

Theorem 1: If W has the non-central Wishart distribution given in (|7j), then the pdf of the latent 
roots A = diag(Ai, . . . , Am.) of |W — AS| = depends only on the latent roots O = diag(a;i, . . . , ujm) 
of |MM^ — ujH\ = 0, and is given by 



m{m—l) 



/a(A) = e-""oFi(n;0,A) e"^-^! Al"- n(^^ " ^.)'- (32) 

Due to the scaled identity matrix structure of the covariance matrix S in our case (see (j3))), it is 
easily seen that the latent root distribution A of the matrix W, as required in (j22j) . can be obtained 
from the distribution given in ()32|) by noting that 2cr^Aj = Aj and 2a'^Ui = uji for i = 1, . . . m, and thus 
2a^A = A and 2(T^r2 = fi, where we have denoted the latent roots matrix of MM^ by fi. Hence, by 
applying this change of variables to we get the required eigenvalue distribution of the matrix W 
as 

/a (A) = (1 + /€)"e-^*%Fi(m; ^A)^^^-^ e-^'^^^I^AI™-" FT ^ _ ^ . (33) 

From the definition of M in (jSJ we see that MM^ = ^i^'^ rn,n^m,n = nfi'^^m where "^rn denotes 
the m X m matrix of all ones. By decomposing as = rnuu^ where u is the m-vector = 
^^[1, !,...,!] and noting that u^u = 1, we observe that the only non-zero eigenvalue of the matrix 

is equal to m. It follows that the only non-zero eigenvalue of the matrix MM^ is equal to mufi'^, 
and thus ui = mujj? and tUj = for z = 2, . . . , m. Substituting these into (jH^ and using the definition 
of the Rician factor k from (jS)) we have 

/,,,...,,„(Ai,...,AJ = (l + «:)-"e-™%Fi(n;(l + /s:)ri,(l + «:)A) X 

^m{m-l) 



-(l+«)E.=iA« 



Trr,{n)Tm{m) 



wk] n(^^-^^)'- (34) 
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Note that, when k = 0, (jMjl reduces to the distribution of the Rayleigh fading latent roots, given in 
[S], since in this case fl = Omxm and oFi{n; Omxm, A) = 1. 

The general capacity expression for the Rician fading channel for an arbitrary number of trans- 
mit/receiver antennas is then given from (j22|) . by. 



CNT,Na = [ y^log ( 1 + J /ai,...,A„(Ai 



Xm)dXl . . . dXra , (35) 



where /ai,...,a™(Ai, . . . , A.^) is given in (jHll) above. 



V. A New Signalling Scheme for Multiple Transmit Antenna Systems in Rician 

Fading 

In Section riV-B.2l we observed that there is a large gap between the capacity upper bound for the 
Rician multiple transmit antenna system under the assumption of known k at the transmitter and the 
capacity of the system without this assumption. In this section we propose a better signalling scheme 
for multiple antenna systems that achieves higher capacity by explicitly making use of the knowledge 
of the Rician factor at the transmitter. 

Recall from Section ITTTl that the optimal signal choice for such a multiple antenna system in Rician 
fading, subjected to an average transmit power constraint P, is zero-mean complex Gaussian. Thus, the 
only thing we do not know is the covariance matrix structure of the optimal input signal distribution. 
Based on the discussion at the beginning of Section 11111 we propose the following choice for the 
covariance matrix of the zero-mean complex Gaussian input signal x: 

Q'^ = .J A lNr+>^^Nr) (36) 

where, as before, "^Nt is the Nt x A^^- matrix of all ones. Note that, when k = 0, (j36|l becomes the 
optimal covariance for the Rayleigh channel. On the other hand, as n — > oo, Q*^ — > Q°° which is 
the optimal covariance for the AWGN MIMO channel. Thus Q"^ reduces to the optimal covariance 
matrices at these two extremes. 

With this choice for the covariance matrix Q, the capacity of the multiple antenna system becomes 

C-Mt,n, = En jlog det (^ ^^^^ ^ H (I^, + k^^,) + I^v.^ } • 

Note that we may write = where e denotes the A^t-— vector of all ones. Since the matrix H 
reduces to an length row vector when N^i = 1, the capacity of the multiple antenna system in this 
case can be written as 



(37) 



where as usual W = HH^ = J2nT=i I (^)i,nr 1^ ^nd we have defined Z = e^H^ = J2nr=i (H)i,„^- 
Since the (H)-^^^'s are independent complex A/'c ^;^(1 + j), 2cr^j random variables, it follows that 

Z ~ Afc (^^r:^(l + j), ^Ntct"^^. Hence, using our earlier notation, it can be easily shown that |Zp ~ 

Wi (1, ^, Ntk) and W ^ (iVr, i^, Ntk) . 

It is clear that these two random variables Z and W are not independent and so we do not know 
their joint distribution, which is required for evaluating fl37|) . Thus we resort to capacity bounds and 
derive both upper and lower bounds on the capacity of the multiple antenna system with the choice 
(j3b|l . In particular, we obtain a tight lower bound on the capacity which shows that the proposed 
choice of the covariance matrix is far superior to the scaled identity covariance matrix for any non-zero 
K (and, of course, is the same as that capacity when k = 0). 
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A. Upper Bound for C^^ ^ 

Applying Jensen's inequality to (jSZD and noting that E{\Z\^} = ^{1 + Ntk) and E{W] = Nt, 
we have the following upper bound on the capacity of a multiple transmit and single receiver antenna 
system in Rician fading with the proposed covariance matrix: 

It can be shown that in the special cases of k = and n — > oo, the upper bound (PH|) reduces 
respectively to, 

< log(l + P), (39) 

and 

C^^^ < hgil + NtP). (40) 

From ()25j) and the results in [3] we know that the right-hand side of ()39|1 is in fact the exact capacity of 
the system in this case as Nt — >■ oo. Hence, in the case of k = the upper bound (jHEI) is achieved as 
Nt — > oo. Also, from the remarks in Section FlV-AI following (I^Hjl . we see that the right-hand side of 
fjlUj) indeed is the exact capacity of the system in this case for any value of Nt- Hence, when k, — > oo 
the upper bound (|HHj) is achieved for any value of Nt- 

B. Lower Bound for C^^ ^ 

Since both W and Z are non-negative random variables we may obtain the following lower bound 
on the capacity of a multiple transmit and single receiver antenna system in Rician fading with the 
proposed input covariance matrix: 

Cn^,i > (! + «:) exp{-NTK) log (yY^z + exp(-(l + K)z)h {2^Ntk{1 + k)z^ dz. 

In Figs. IHlandiniwe have shown the derived bounds for the capacity of a multiple-transmit antenna 
system along with the capacity corresponding to the scaled identity covariance matrix. From the lower 
bounds shown on these figures it is clear that the proposed signalling scheme achieves much higher 
capacity than that of the scaled identity covariance matrix for sufficiently large values of n or Nt- It 
is also observed that the difference between the upper and lower bounds decreases as Nt increases. 
From Fig. El we note that when /t is large the upper and lower bounds are almost the same unless 
the number of transmit antennas is very small. Thus, for large k, a reasonable approximation to the 
capacity of the proposed scheme can be obtained by taking the large-K asymptote of the upper bound 

ck,i ^ ^^^{^ + ^T$) ^^^^ 

In Figs. IHl and IHl we have also shown this large k, approximation to the upper bound of the capacity. 
We observe that indeed (PT|) is a very good approximation to the capacity even for relatively small 
values of k- Note also that (^1]) gives the exact capacity in the cases of k = and k — > oo. 

Finally, it is worth noting that in these figures we have also included the capacity upper bound for 
a Rician channel with receiver CSI derived in Section UTTl From Fig. |Hlwe observe that for relatively 
small values of k there is still a significant gap between the general upper bound for the Rician channel 
in this case given by (jSDl) and the upper bound on the capacity of the proposed new design given 
by (j38j) . However, as k increases we observe from Fig. El that this difference also becomes smaller, 
although the proposed scheme still does not achieve the upper bound (j2(J|) . 
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VI. Conclusions 

In this paper we have investigated the average-power-constrained capacity of multiple antenna sys- 
tems under Rician distributed fading when the receiver has access to channel state information, but 
not the transmitter. We have considered two different cases concerning the knowledge of the fading 
available at the transmitter: that in which the transmitter has no knowledge of the fading at all; and 
that in which the transmitter has knowledge of the Rician factor k but not the exact value of CSI. 
While obtaining the exact capacity in the former case, we were able to derive lower and upper bounds 
for the latter case. The exact capacity in the former case also quantifies the capacity variation of a 
multiple antenna system designed to be optimal for a Rayleigh fading channel but in fact operating 
in a Rician fading environment. For this case, we derived an integral expression for the capacity of a 
general system having an arbitrary number of transmit /receive antennas. In some special cases, we nu- 
merically evaluated this capacity expression. We specifically investigated the capacity of such systems 
for large numbers of transmitter antennas and when only one end of the system (either transmitter or 
the receiver) is equipped with a multiple antenna array. 

A new signalling scheme has been proposed for the case when the transmitter knows the value of 
K, though not the exact CSI. We have analyzed the capacity of this new scheme, in terms of lower 
and upper bounds, for a multiple transmit antenna system, and demonstrated that it offers much 
higher capacity than that of the unknown-K capacity-achieving distribution. We also derived a simple 
approximation to the capacity of this scheme for sufficiently large values of k. 

Our results indicate that Rician fading can improve the capacity of a multiple antenna system, 
especially if the transmitter knows the value of k. Moreover, the proposed signalling scheme provides 
a mechanism for capturing this improvement. 

Appendix 

In this appendix we present a few mathematical concepts that have been used throughout this paper. 
Most of these are related to various types of special functions needed in our analysis. 
The complex multivariate gamma function, Tm{a)i is defined as |7j 

„ m 

fm{a) = / e-*'^[^l|A|'^-'"(rfA) = 7r^'"('"-i)TTr(a- (A;- 1)). (42) 

Jah=a>o t=i 

Note that, it follows from ^ that f i(a) = T{a). 

The generalized hypergeometric function ^2), [Zli is defined as 

pFg{a^,...,ap;h,...,bg;x) = 2^jj-. TT^Ti (43) 

^ {Ol)k ■ ■ ■ [Oq)k kl 

where p and q are integers, and the hypergeometric coefficient {a)k is defined as the product 

(a)fc = a{a + l)...{a + k-l), (44) 

with (a)o = 1. 

The complex hypergeometric function pFg{ai, . . . ,ap;bi, . . . ,bq; S) of an n x n Hermitian matrix S 
can be defined as [Zj 

(ai)^ . . . (ap)^ C^{S) 
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where p and q are integers, Ck{S) is the zonal polynomial (JD], ^1], ^H], JE]) of the n x n Hermitian 
matrix S corresponding to the partition k, = {ki, k2, . . . , kn), ki > k2 > ■ ■ ■ > kn > 0, of the integer k 
into not more than n parts and [0]^ is the complex multivariate hypergeometric coefficient defined as 

n 
i=l 

Note that Ck{S) is a homogeneous symmetric polynomial of degree k in the latent roots of the matrix 
S. 

The hypergeometric functions with two argument matrices, S and T (both n x n), can then be 
defined via the relation [H], [Z|, [IS], 

pFq{ai,...,ap;bi,...,bq;S,T) = / .pFq{ai, . . . ,ap;bi, . . . ,bq; SUTU^){dU) 

JU(n) 

where U{n) is the unitary group of all ri x n complex unitary matrices U , i.e. UU^ = In and (dU) is 
the invariant (Haar) measure on U{n) normalized to make the total measure unity. 

A special case that we need is 0-^1 , which is the Bessel function of matrix argument ^Zj, ^H], |Zj, 
which can also be represented as 

oFi(m;M^) = [ e''^'''''+^'\dU) (45) 

Ju{m) 

where if is an n x m complex matrix with n < m, U = [Ui, U2] with Ui being an n x m complex 
matrix and Uf^Ui = In- Note that, in H and (dU) denote the complex conjugate of the matrix 
H and the normalized invariant measure on the unitary group U{n), respectively. 
It can be shown that for a scalar r, (j45|) reduces to 

oFi (m;r2) = r(m)r-("^-^)/„„i(2r) (46) 

where Iu{z) is the u-th order modified Bessel function of the first kind jS]. 
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Asymptotic Capacity for Large Number of Transmit Antennas Versus Rice Factor 
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Fig. 1 

Asymptotic Capacity for Large Number of Transmit Antennas Versus Rician Factor. SNR = dB. 
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Asymptotic Capacity for Large Number of Transmit Antennas Versus Rice Factor 
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Fig. 2 

Asymptotic Capacity for Large Number of Transmit Antennas Versus Rician Factor. SNR = 10 dB. 
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Asymptotic Capacity for Large Number of Transmit Antennas Versus Rec. Antennas 
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Fig. 3 

Asymptotic Capacity for Large Number of Transmit Antennas Versus Number of Receiver 

Antennas, k = 10. 
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Rician Channel Capacity for Single Receive Antenna Versus Transmit Antennas 
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Fig. 6 

Rician Channel Capacity for Single Receiver Antenna Versus Transmit Antennas. 
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Rician Channel Capacity for Single Receive Antenna Versus Transmit Antennas 
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Fig. 7 

Rician Channel Capacity for Single Receiver Antenna Versus Transmit Antennas. 
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Capacity of the New Scheme for a Multiple-Transmit and Single Receiver Antenna System 




CO 5 



ra 4 
O 



cc 

? 3 



as 

to 2 
O 



-*- Rayleigh Signals 

Capacity Upper Bound 
New Capacity Upper Bound 
New Capacity Lower Bound 
Upper Bound Aymptotic 



5 10 
Number of Transmit Antennas 



Fig. 8 

Capacity of a Multiple-Transmit Antenna System in Rician Fading with Proposed New Signalling 

Scheme. Nr = 1 and k = 1. 
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Fig. 9 

Capacity of a Multiple- Transmit Antenna System in Rician Fading with Proposed New Signalling 

Scheme. A^^ = 1 and k = 10. 
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